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Abstract 

We investigate the performance of a large class of cosmic microwave background experiments with 
respect to their ability to measure various cosmological parameters. We pay special attention to the mea- 
surement of the total cosmological density, Q. We consider interferometer experiments, all-sky single-dish 
experiments, and also single-dish experiments with a deep-patch technique. Power spectrum estimates for 
these experiments are studied, and their induced errors in cosmological parameter estimates evaluated. 
Given this motivation we find various promising corners in the experiment parameter space surveyed. 
Low noise all-sky satellite experiments are the expensive option, but they are best suited for dealing with 
large sets of cosmological parameters. At intermediate noises we find a useful corner in high-resolution 
deep patch single-dish experiments. Interferometers are limited by sample variance, but provide the best 
estimates based on the very small angular scales. For all these experiments we present conservative, 
but still promising estimates of the accuracy of the measurement of f2. In these estimates we consider a 
variety of possible signals not necessarily in the vicinity of standard cold dark matter. 



1 Introduction 



In many recent papers |l], ||, ||, [|, || it has been shown how the cosmic microwave background (CMBR) 
angular power spectrum Ce may provide a clean measurement of several cosmological parameters, and also 
of the so-called inflationary observables. In particular the Universe total density in units of the critical 
density fl, is known to leave a distinct imprint in the Ct spectrum ||. Therefore it has been suggested j|, ^ 
that the next generation of experiments, in particular MAP [Q and COBRAS/SAMBA ||, should settle 
the question of the total density of the Universe. This would also determine the geometry of the Universe, 
deciding at long last which of the three types of Freedman models (spherical, planar, hyperbolic) describes 
our Universe. A standard error bar in f2 of about 10% is expected for an experiment with features similar 
to what is now known as the MAP proposal |Q. In this estimate one marginalizes with respect to a large 
parameter space, with unassuming priors. 

In this paper we reexamine and extend the work in || where estimates for cr(f2) were computed for 
prototype satellite experiments. We briefly assess the errors and uncertainties in the covariance matrix 
approach used in |Q, |). We find that the results obtained by this method are to be seen as a guide only. 
Slight changes in the algorithm, as well as in the number of parameters considered, may easily raise a 
standard error in fl from, say 10%, to about 30%. Also in the analysis performed in ||, ||] it was assumed 
that the "signal" comes from a standard cold dark matter theory (sCDM). We find the unsurprising result 
that if one considers instead the case of signals coming from theories with a scalar tilt smaller than 1, or low 
f2 models, then the errors are significantly larger. Since we don't a priori know what the signal is, one must 
be prepared for the worse when estimating the performance of a given experiment. Overall we find that the 
analysis in j|, ^ is a perfectly valid preliminary estimate. However, when reinserted in the context of its 
own uncertainties, and reapplied to less generous signals, this analysis reveals that the conclusions in |4|, |^| 
are perhaps embedded in what might be described as realistic optimism. In this paper we take the opposite 
point of view: we decide to be pessimistic, not to be cumbersome, but instead in order to investigate how 
one could vary the observational strategy in order to improve final estimates of f2. 

We consider a covariance matrix set up yielding the worst prediction for sCDM, and also consider the 
case of signals coming from non-sCDM theories which make more difficult to measure. However we also 
generalize considerably the space of possible experiments covered by the analysis in 0, ^. Taking up work 
in |], [To] we consider single-dish experiments in which the sky-coverage is deliberately small (the so-called 
deep patch technique). We also consider the case of interferometric experiments (see Jl3| for a good review). 
We then try to design the ideal experiment for measuring f2, given a constraint mathematically rephrasing 
fixed finite funding. This work provides guidance along two lines. Firstly, one may examine what is the 
ideal scanning method for best results subject to the constraint. Secondly, one may provide the value of the 
error in f2 as a function of this constraint, assuming ideal scanning. This will set up lower bounds on the 
constraint for a meaningful experiment, telling us also how fast these errors will go down thereon, after a 
given constraint improvement. Not surprisingly we find more than one promising corner in the experiment 
parameter space surveyed. In summary, all-sky low noise satellite experiments, high resolution intermediate 
noise ground based experiments, and interferometers, all, for different reasons, provide estimates for f2 with 
useful errorbars. It is now up to political details to decide who will win the race of the weighing of the 
Universe. 

Our paper is organized as follows. In Section [2] we review the covariance matrix approach to the estimation 
of cosmological parameters. However we tie this approach to previous work we have done [jlo| on power 
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spectrum estimates for a large class of experiments. Then in Section |j| we make use of this formalism to 
produce heuristic arguments on what are the best scales for measuring various cosmological parameters, with 
special reference to Q. We then show how the scanning strategy may improve experimental performance 
on sections of the spectrum which are most relevant for the determination of tt. We show cases where 
reducing sky-coverage is useful for improving the measurement of fl, and cases where it is not. Given these 
insights we then jump into the full problem of estimating the errorbar in £1 given the specifications of an 
experiment ideally suited to measure subject to its constraints. In Section |] we present results for signals 
coming from the vicinity of sCDM, and in Section [7] we consider other signals. In Section [| we also spell 
out formalism uncertainties and how they affect the results in Section |[ As a side remark in Section || 
we investigate the interferometers' performance with respect to other cosmological parameters as well. We 
conclude with a few comments on what steps should be taken towards a more precise and comprehensive 
analysis of determination of cosmological parameters. 

2 Determination of cosmological parameters with CMB experi- 
ments 

We start by reviewing the covariance matrix approach j^, 11, |l2) in the guise to be used in this paper. To 
set the notation, the Cg spectrum is defined from the two-point correlation function C{9) as 

C{6) =]T^-QP,(cos0), (1) 
or alternatively from Ci = (|a^ m | 2 ), where the ag m are the spherical harmonic coefficients in the expansion 

^P = f; J2 ( 2 ) 

1=0 m=-t 

The Ci spectra will here be used to estimate a maximal parameter space spanned by s = {Q, f2, flth 2 , h, ns, A, 
r, tit, t, Here O, fif,, and f2„ are the total, the baryon, and the heavy neutrino densities, in units of the 
critical density. When fl v ^ we have assumed one massive neutrino, h parameterizes the Hubble constant 
in the usual way: H — lOOft-Km sec _1 Mpc _1 . A is the cosmological constant. Whenever A ^ we have 
defined the total density of the Universe, responsible for its geometry, to be ft = SIq + A, where flo is the 
total density in matter, radiation, and neutrinos. n$ and tit are the tilts in the primordial power spectrum 
of scalar and tensor perturbations, respectively. Q is the quadrupole in [iK , so that Ci = (An/5)(Q /T ) 2 , 
where To is the average CMBR temperature. We shall normalize all spectra by the 10° rms temperature, r 
is the ratio of scalar and tensor components defined at £ = 2, that is r = Cj/Cf • Finally r is the optical 
depth to the epoch of recombination. Besides this maximal set of parameters we will also consider subsets 
in order to emphasise the dependence of the final results on the number of uncertain parameters. 

We shall use a Gaussian approximation to the likelihood function as proposed in ||, [ll], |l2) . This consists 
of approximating the likelihood linearly near its maximum, resulting in 
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where Sq is the underlying theory set of parameters, and the curvature matrix [a] is given by 

-dCtMdCifo) 



Y — 



dsi dsj 



(4) 



Inverting [a] one therefore obtains the parameters covariance matrix [C] = [a]^ 1 , thereby converting the 
errors in the measurement of the Ci (given by <r 2 (CV)) into errors in the estimates of s. The diagonal 
components, Ca, in particular, are the errors er 2 (si) in the parameter s, obtained by marginalizing with 
respect to all other parameters. We will find that the marginal errors a 2 (si) depend very strongly on the 
total set of parameters chosen. 

The errors in the estimates of the power spectrum cr(CV) are due to cosmic/sample variance, instrumental 
noise, and foreground subtraction. Following [|l3| [l4|] we shall translate the effects of foreground subtraction 
into a renormalized noise term, that is we take foreground subtraction errors into account by multiplying 
the instrumental noise by a deterioration factor. The deterioration factor is a priori unknown. An optimist 
would place it between 1 and 2, a pessimist above 3. We have empirically found this procedure to be 
appropriate for interferometer experiments Jig] by directly examining the effects of deconvolving galactic 
foregrounds in real data using maximum likelihood methods. We shall consider errors in Ci estimation for 
a large class of experiments. These include all-sky single-dish experiments, single-dish experiments with 
a deep-patch technique, and interferometers. Power spectrum estimate errors for these experiments were 
studied in ]9|, [l^] , where the reader may find more detail and derivations. Here we merely review the main 
results in g, |o|. 

2.1 A satellite experiment 

We first review the case of a single-dish experiment covering a large portion of the sky / s k y (see H |[ [l(J for 
more details). Let the pixel area be f2 p ; x = 7r# 2 , where 6^ w 2.35crf, is the FWHM size of each pixel, assumed to 



be Gaussian shaped. The probed spectrum will then appear multiplied by a beam factor B(l) — exp(P/2of 



which needs to be deconvolved. Now let the noise per pixel be i7p ix = s /i p i x , where s is the detector 
sensitivity and i p i x is the time spent observing each pixel. Then fixing the detector sensitivity s and total 
time of observation i t ot fixes the quantity u> -1 = 47rs 2 /t to t = fpix^pix/ fsky, where / s k y is the fraction of sky 
coverage. We therefore parameterize the noise level with w . In this argument we have assumed active 
incomplete sky-coverage, that is, an experimental strategy in which one deliberately focuses on an incomplete 
patch of the sky. Because this is done deliberately, the integration time on each pixel is larger. Therefore the 
errors due to noise are smaller, even though the errors due to sample variance are larger. This set up is also 
sometimes called deep-patch technique. On the other hand, passive incomplete sky-coverage occurs if one 
observes the whole sky, but then portions of it have to be thrown out for various reasons, eg. point source 
contamination, or galactic obscuration. In such a case the noise parameter is w^ 1 — 47rs 2 /t to t — ""pix^pix- 
Then using the results in |g one may prove that 
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for active incomplete sky-coverage. For passive incomplete sky-coverage the appropriate formula is 
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2.2 A deep-patch single-dish experiment 

These formulae break down if / s k y <C 1. The effects of small sky coverage have been studied in H [h]] and can 
be summarized as follows. Firstly the probed spectrum becomes the convolution of the raw spectrum with 
the window spectrum. For definiteness let us assume that the field is a square with size L, in radians. Then 
the probed spectrum will differ from the raw spectrum in that features on a scale Al w 1/L are smoothed 
out, and also a white noise tail appears for I < 2tt/L. If the distortions introduced are considerable, then a 
deconvolution recipe becomes necessary, increasing the errors significantly. We have however checked that 
sCDM Doppler peaks do not require deconvolution for any window with L > 4° , as long as a bell is applied 
to the field, if this has edges. The I < 2tt/L section of the spectrum cannot be recovered, but apart from this 
the probed and the raw spectra are proportional to each other. If Q < 1 the lower bound on L for avoiding 
the need for deconvolution becomes even less restrictive (of order L > 2° for fl = 0.3). 

Secondly incomplete sky coverage introduces correlations between the modes used in the Ci estimates. 
This reduces spectral resolution, that is, it allows independent Ci estimates only with a separation Al > 1 
(see JlCj for formulae for Al). Correlations also reduce the number of independent modes available for each 
independent estimate, and therefore cosmic variance is increased by finite sky-coverage (an effect sometimes 
called sample variance). These two effects can be studied quantitatively by stereographically projecting the 
observation field onto a plane, Fourier transforming, and replacing the Fourier plane by what we called an 
uncorrelated mesh of modes. This mesh contains only quasi-uncorrelated modes, and their density equals 
the density of independent modes present in the full Fourier plane itself. For a square field with a side 
L treated with a cosine bell the side of the uncorrelated mesh is fco ~ 2-k/L. From the uncorrelated 
mesh modes k; one may then extract power spectrum estimates for values of U for which mesh points exist 
satisfying U — int(|k;|), where int denotes the integer part. In general one may come up with uncorrelated 
estimates only with a separation Al w ^/l 2 + k^/it — I, for I > fco- Only for I > /cq/(27t) can individual Ci 
be estimated (Al ss 1). The curvature matrix [a] corresponding to these power spectrum estimates is then 
given by 



where N(li) is the number of modes satisfying k = int ( | | ) . 

If there is near all-sky coverage then At ss 1 and N(l) = (21 + l)/ s ky, recovering the large sky-coverage 
formulae. In fact the uncorrelated mesh results for L w 202° (corresponding to / s k y ~ 1) differ from the large 
sky-coverage limit results by less than 1%. On the other hand for L < 100° (that is for / s k y < 0.25) the errors 
inferred from the large sky-coverage limit are always grossly underestimated. Therefore the uncorrelated 
mesh formalism acts as a generalization of the large sky-coverage formulae. 

2.3 Interferometers 

For interferometers the field is Gaussian shaped (the so-called primary beam [|l3| Q). Let d w = 2.35a w 
be the primary beam FWHM. Deconvolution problems may be avoided for sCDM by imposing a primary 
beam size 6 W > 4° (9 W > 2° for f2 = 0.3). Interferometers make measurements directly in Fourier space 
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where the errors in the estimates a 2 (Cn i ) are 
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(which in the community jargon is called the uu-plane). Again uncorrelated estimates and their variances 
may be obtained by means of an uncorrelated-mesh |L0|, this time with a side ko » 2^/Att \og2/9 w . Also if 
one decides to observe n/ well-separated fields, then each mesh-point acquires an extra index i = 1,. . . 
and points with different indices are uncorrelated. With these modifications the curvature matrix may then 
be computed as in (^). 

Now let A v ; s be the number of visibilities in each uncorrelated mesh cell, let s be the sensitivity of 
the detectors, and t V i S be the time spent observing each visibility. The coverage density is given by p c = 
N v i s t V i S Q, w /tt, where tt is the time spent on each field and Ct w — ira^ is the field area. We shall assume 
that the coverage density is uniform in a ring of the uw-plane going from k to a certain k max delimiting the 
outermost uv-tracks depicted by the interferometer. This may be attained with a dish geometry like the one 
proposed in . 

In order to parameterize the noise for interferometers, we now notice that for fixed detector sensitivity 
and total observation time one should now keep constant u> -1 = (2ir) 2 s 2 / (p c ttot) — (^ 7r ) 2(J N / (^ s3n f)- This 
allows comparing experimental strategies subject to the same constraint (time and money) but choosing, 
say, different primary beam sizes and number of fields. As shown in Jic[| the power spectrum estimates for 
interferometers provided by their uncorrelated mesh are now affected by the errors: 



It should be remarked that as we go up in I the signal-to noise in each independent mode goes down as a 
power law (we are assuming that the outer ring is before the Silk damping tail). This is to be contrasted 
with single-dish experiments, where the same signal-to-noise always goes down like an exponential in I. This 
feature makes interferometers desirable for the measurement of high I features. 

3 Relevant scales for determining cosmological parameters 

We now investigate which are the statistically most relevant Ce's for the determination of the various cos- 
mological parameters. To do this we start by assuming for each parameter Si that all other parameters are 
kept fixed, so that 

1 ^S(l; Si ) 



\^ 2*1*1 (io) 



a 2 {si) ^ N{£) > 



defining a sensitivity factor S(l; Si) 



S(l;si) = 



dC e (s Q ) 
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and a noise factor N{€) which for a single-dish experiment with active limited sky coverage is given by 



m) = / s - ky 



C t B{lf 



(12) 



The sensitivity factor S(l; Sj) tells us the relevance of each Cg for the determination of the parameter Si 
subject only to the unremovable cosmic variance. The noise factor N(£) will then deteriorate this result. 
Knowledge of where S(l; Sj) is largest will then suggest experimental features. If one is most interested in a 
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Figure 1: The Ci spectrum of sCDM (line) and of theories in its vicinity obtained by changing one of the 
parameters Sj to a value indicated in the label (points). The changes introduced may be mild or radical 
depending on the parameter. Also there is some level of degeneracy in the changes introduced. 
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Figure 2: The sensitivity functions S(l; Si) for various cosmological parameters as inferred from the variations 
used in Fig.l. We have taken percentage sensitivities S(l; Sj)s? on the plots on the left hand column. 

given parameter then one should require that N(£) be the closest to 1 (its mathematical minimum) where 
S(l; Si) is the largest. In practice the determination of a given Si is coupled to the determination of all the 
other parameters. Still one may detect general trends by studying the S(l; Sj) curves. 

3.1 The sensitivity to cosmological parameters 

The Ct spectrum dependence on s is illustrated in Figure [j]. We have assumed that we are in the vicinity 
of sCDM defined as s = (18,1,0.01,0.5,1,0,0,0,0,0). We have then shown spectra obtained by changing 
various parameters s% keeping all other parameters fixed at their sCDM value. 

Two comments should be made. Firstly the function Ce (s) changes at rather different rates depending on 
the parameter s, . The scalar tilt ns and the total cosmological density f2 change the whole spectrum rather 
drastically when subject to small percentage variations. A multiplicative shift in I is the benchmark of Q. 
An overall tilt of the spectrum is of course the hallmark of ns- Late reionization (r ^ 0) and tensor modes 
(r ^ 0) affect strongly the relative normalizations of the plateau and peaks. The separate shapes of these two 
spectrum sections, however, remain roughly unchanged. Hence a low coverage high resolution experiment, 
and a high coverage low resolution experiment, considered on their own would effectively be blind to r and 
r. The other parameters considered leave rather subtle imprints. The baryon content enhances odd peaks 
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relative to even peaks, the so-called "acoustic signature". A large Hubble constant reduces all peak heights. 

As a second remark we note that changes in different parameters may change certain sections of the 
spectrum in the same way. In particular if one only has access to the primary peak a very high level of 
degeneracy will occur. This degeneracy will then correlate estimates of degenerate parameters, increasing 
their marginal variances significantly. Again fl and n$ are the only two parameters which appear not to be 
much affected by this problem, as their signatures are sufficiently distinct. For disentangling h and Vl^h 2 , 
on the other hand, one should go into the secondary peak region. In particular the second Doppler peak is 
essential for removing the degeneracy between these two parameters. For all other parameters one has to 
rely on rather subtle differences in the spectrum shape to remove their degeneracy. 

These statements are formalized in Fig. ||where we have plotted percentage sensitivities (that is S(l; Sj)sf ) 
wherever Sj ^ 0, plain sensitivities otherwise. Clearly most of the information leading to f2 and the fixing of 
other parameters which may have a degenerate effect (like h) is contained in the Doppler peaks, particularly 
in the secondary Doppler peaks. The plateau carries little relative weight. The sensitivity S(l;Cl) has an 
oscillatory structure since most of the effects of changing SI consist of a shift in I. Hence around stationary 
points of the spectrum there is no information on Q whereas the information reaches peaks near the steepest 
slopes of the spectrum. Since the shift is multiplicative, the change in the spectrum is largest at high I. 
Hence the secondary Doppler peaks' slopes, rather than the main peak summit, are the ideal scales for 
measuring fl. For h and flbh 2 most of the information is in the summit values of the Doppler peaks. The 
sensitivities to ns, r, and r are smoothly distributed all over the spectrum. Overall, the sensitivities for il 
and ns are much larger than for other parameters. 

From Fig. |^ one may also draw information on degeneracy. This may be inferred from similarly shaped 
S(l; Si) curves. In the absence of noise: 

^ M = ^^^^-. d3) 

Hence sections of the spectrum where S(l; Si) oc S(l; Sj) contribute to a large correlation between the Si and 
Sj estimates. The fact that S(l; SI) looks so different from all other curves then gives credit to the hope that 
a clean unambiguous measurement of SI may be achieved by a CMB experiment. 

3.2 Noise factors and scanning strategies 

The results of the previous subsection open the problem of the ideal scanning strategy for measuring cos- 
mological parameters. Suppose that we have a single-dish experiment with active finite sky-coverage. For a 
fixed resolution and noise constraint w~~ l one now has to choose the best coverage area. Best results involve 
a delicate balance between noise and sample variance. Small fields have an increased sample variance, and 
are blind to the low I section of the spectrum. However they enable reducing the effects of noise signifi- 
cantly, allowing for more integration time to be concentrated on each pixel. A better measurement of high 
I sections may be then provided. The exact balance point between these opposite requirements depends on 
the particular question one wants to answer. For the determination of cosmological parameters this balance 
point is set by demanding that N{t) be small where S(l,Si) is large. In the previous subsection we showed 
which scales I should then be targeted. 

In Figure [| we show the effect of reducing the sky-coverage in two experiments: one with low noise 
and intermediate resolution, the other with intermediate noise and high resolution. As the sky coverage is 
reduced, the low I section of the spectrum becomes more noisy, but after a certain cross-over / one actually 
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Figure 3: The noise factor N(£) for sky coverage from L = 200° to L = 20°, for two experimental constraints. 
On the left w" 1 = (15^iT) 2 (deg) 2 and 6 b = 0.5°. On the right w~ l = {60fiK) 2 (deg) 2 and 9 b = 0.2°. 

obtains better estimates. In the first experiment the cross over occurs at I as 400, when the noise factor has 
already shot up exponentially. Therefore in this case it does not seem that reducing sky coverage may be of 
any use. The real limitation is resolution and not noise. One is better off mapping as much sky as possible 
in order to reduce sampling effects in the I < 300 sections of the spectrum, where the experiment fares well. 
In the second experiment, on the contrary, by reducing the sky coverage we throw away the I < 200 section 
of the spectrum, but we also greatly improve results at high I. Reducing sky-coverage in this experiment 
allows for noise to be defeated effectively, letting the experiment high resolution manifest itself in useful 
estimates of I > 400. Due to the high sensitivities in in this section of the spectrum one may then achieve 
an accurate measurement of fi. The results from these two experiments may then be of comparable quality. 
However this quality is achieved by targeting different sections of the spectrum, and with rather different 
ideal scanning strategies. 

These two situations are typical. In the first case one relies on a very accurate mapping of the first peak, 
by means of low noise, large sky coverage, and intermediate resolution. In the second, on the contrary, one 
makes use of the large sensitivity f2 has at large I, and targets this section only. Low coverage, and high 
resolution are the right combination to do this at intermediate noise. 

Although one can never give up the uncorrelated mesh approach in the case of interferometers, their 
noise factor may be approximated by 



In Figure ^ we have plotted N(i) for an interferometer with a primary beam of 9 W = 4°, and w^ 1 = 
(10/iisT) 2 (rad)~ 6 (left) and w^ 1 = (30/iif ) 2 (rad)~ 6 (right). The various curves correspond to different 
numbers of fields. The most noticeable feature is that although the overall noises are high, they do not 
increase by much with £. As a result interferometers provide the lowest N(£) at high £. The main limitation 
is clearly sample variance, especially if the noise is of order w^ 1 = (10^i4T) 2 (rad)~ 6 . For this noise level the 
best scanning strategy will always be to try and map as many fields as possible, granted the limitations of 
parallel point source subtraction. Sample variance becomes less relevant at very high £ however, and so even 
with a limited number of fields one may fare well at high £ with interferometers. Once again, given the high 
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Figure 4: N(£) for an interferometer with a primary beam of 9 W = 4°, and w^ 1 = (10/iif ) 2 (rad)~ 6 (left) 
and u> -1 = (30/iif ) 2 (rad) -6 (right). The various curves correspond to different numbers of fields starting at 
1 (top curves) up to 19, in steps of two. 

sensitivity of the spectrum to Q at high I one may expect a good measurement of fl with an interferometer. 
Hence interferometers may be expected to be competitive with the previous two experimental strategies, for 
yet another reason, and with yet another ideal scanning strategy. 

4 Ideal CMB experiments for weighing the Universe near sCDM 

The heuristic arguments of the previous Section will now be replaced by a full calculation of errors in 
for a large section of experiment parameter space. As stated before best results in any CMB experiment 
involve a delicate balance between noise and sample variance. The exact balance point depends on the 
particular question one wants to answer. In || we considered the problem of detecting secondary Doppler 
peaks. In comparison, we have found that the problem of the measurement of f2 marginalizing with respect 
to a large set of paramaters is more sensitive to sample variance, and therefore always requires a larger 
optimal coverage area. If the set of uncertain parameters is reduced, however, not only do errors go down 
sharply, but also sample variance again becomes of secondary importance. We will illustrate this point in 
Section ^|. Unless otherwise stated, we will, in this section, marginalise over all the parameters s = {Q, 57, 
flbh 2 , h, ns, A, r, ut, t, 0^]-, except for fl u which we assume to be zero, i.e. we assume no massive neutrinos. 

We will highlight three corners of experiment parameter space which we found the most promising: 
satellite experiments, ground based deep-patch experiments, and ground based interferometers. 

4.1 Satellite experiments 

An instrumental noise level of w^ 1 = (7.5/iK) 2 (deg) 2 is often quoted as realistic for a satellite experiment 
(H H U)- Some of the COBRAS-SAMBA channels achieve a factor of 5 improvement on this m. Whatever 
the case one must bear in mind that these instrumental noises are then considerably deteriorated by the 
need to deconvolve foregrounds. For definiteness we shall here consider the following two noise levels as 
renormalized after foreground subtraction: w^ 1 = (15/iAT) 2 (deg) 2 and w^ 1 = (30/iir) 2 (deg) 2 . We believe 
that either figure is realistic, but the first corresponds to an optimistic hope, the latter to a pessimistic 
expectation. A remark on normalization should be made. If in future experiments Q comes out very 
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Figure 5: Percentage standard errors in £1 for w 1 = (15fiK) 2 (deg) 2 (top) and w 1 = (30^X) 2 (deg) 2 
(bottom). In both cases, isolines are plotted from 5% to 100% at 5% increments. 

different from COBE's Q = 18fj,K, then this has the effect of renormalizing w^ 1 by the same factor. Hence 
a lower normalization would boost u> -1 unpleasantly. 

We have plotted the percentage errors in O for these two cases in Fig. ||. Clearly, even considering active 
limited sky coverage, the ideal coverage area is all-sky for the u> _1 = (15/iA") 2 (deg) 2 case, and also for 
w^ 1 = (30fiK) 2 (deg) 2 , except perhaps for very small beamwidths 9b < 0.1°. Measuring as one of a large 
set of unknown parameters is obviously very sensitive to sample- variance. Even for noise levels that are not 
too low, like the ones considered, we are dominated by the need to reduce sample variance, not noise. 

As explained further in Section || there are ambiguities in the covariance matrix approach which allow 
for some uncertainty in the errors. Thus, in the calculations presented here, we present the most pessimistic 
results. This resulting error in f2 is of the order of 20% for w^ 1 = (15^iX) 2 (deg) 2 and 9b = 0.5°. These 
errors are quickly reduced to about 5 — 10% with an improved resolution of 9t, = 0.2°. Therefore a satellite 
experiment is always ideally scanned with the largest possible sky coverage, and the final result benefits 
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Figure 6: Percentage standard errors in SI for w^ 1 — (60fiK) 2 (deg) 2 . Isolines are plotted from 10% to 100% 
in 5% increments. 

enormously from improving the resolution. 

4.2 An intermediate noise deep-patch experiment 

If one is observing from the ground, the noise level is always larger. This is because the atmosphere, as 
well as ground spill-over, acts so as to deteriorate the experiments sensitivity often by a large factor. If, 
however, one takes advantage of the fact that a large aperture is no longer a problem, then one may make 
up for this with increased resolution. The ideal coverage area is now very far from all-sky, but, in any case, 
a small patch of the sky is more feasible from the ground. We illustrate this promising corner of experiment 
parameter space in Fig. ^|. We have assumed w~ x = (60/iK) 2 (deg) 2 . A error in SI of around 10 percent may 
then be achieved with a resolution around 0.1°, with a coverage area corresponding to L = 20 — 40°. It is 
clear from the figure that, for such an experiment, increasing the coverage area is extremely undesirable. 

4.3 An intermediate noise interferometer experiment 

In Fig. ^ we show the percentage error in SI for an interferometer with a primary beam of 9 W = 4°. We have 
assumed that a uniform uu-plane coverage up to £ = 1000 has been achieved by means of a sufficient number 
of elements, and their arrangement in a suitable geometry. With these assumptions we have then allowed 
the noise parameter as defined for interferometers to vary between zero and (100/iif ) 2 (rad) -6 . Finite 
sky-coverage errors due to small size and limited number of fields are now substantial. The sample variance 
in the estimates of SI for the interferometer considered can be read off from the w; -1 = axis. Although 
these errors are not small, their limiting values may be approached even for reasonably large values of w . 
The need for point-source subtraction typically limits the number of fields, so unfortunately one may not be 
able to obtain the best result for a given level of noise 

We would expect the next generation of CMB interferometer experiments, such as the VSA, VCA and 
CBI to attain noise level of below w" 1 « (10/i-fQ 2 (rad)~ 6 , for the ww-coverage assumed above, as compared 
to ui" 1 « (200/iX) 2 (rad) -6 for the CAT experiment, after foreground subtraction. Assuming a fixed obser- 
vation time of about 2 years, one could hope to determine SI to within about 15-20 percent by observing 
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roughly 20 fields. 



5 Dependence on number of undetermined parameters 

As one might expected, the number of unknown parameters in a given model has a profound effect on 
the overall error in determining O (or any other parameter). In this Section we illustrate this important 
effect, and also highlight some weaknesses in this method of analysis used, particularly in relation to the 
computation of the curvature matrix. 

5.1 Margins of error in covariance matrix estimates 

The covariance matrix formalism suffers from an ambiguity in the definition of the derivatives dCi/dsi. 
These have to be computed by means of finite differences, which can never realistically be very small. A 
particularly difficult derivative is dCi/dVL. This may be inferred from the approximate formula, valid in the 
vicinity of sCDM 

c(t, n, h, n b h 2 ) « c(en 1/2 , 1, mi 1/2 , n b h 2 ), (15) 



resulting in the relation 



dCj 



^ida hdCh 
n=i ~ 2 d£ + 2 dh ( ' 



The effects of 51 on the plateau can be fitted by a multiplicative factor || which has zero derivative at fl = 1, 
that is, it is a saddle point. Hence the low I signatures of low f2 models do not appear at all in the linear 
approximation, in the vicinity of sCDM. Making use of the prescription ( [l6| ) always results in larger errors 
than say simply computing the derivative with AQ — .05. This is not at all surprising, as even such a small 
difference in £1 results in a difference in I of the order A£ w I Ail/2. This corresponds to taking a finite 
difference Al in ([l6]) which varies with I and is always much larger than 1 at large £. In general increasing 
the size of the differences used in the computation of the derivatives has the effect of reducing the errors. 
The reason we bring up this numerical problem is that it may have a physical implication. If the errors are of 



13 




Figure 8: Percentage standard errors in O for an all-sky experiment with w^ 1 = (15/j.K ) 2 (deg) 2 , as one adds 
more parameters into the problem. The lowest curves consider only (Q, Q) to which parameters are then 
added in the same order as in the full vector s = {Q, Q,, Slfcft 2 , h, n$, A, r, tit, t, 

the order of 10% then it may make more sense to use large finite differences, rather than carefully computed 
derivatives. Therefore the results presented in the previous sections (where we used very fine grids and also 
(|l6|)) are to be seen as conservative estimates. A non- linear analysis will always provide smaller errors. 

In Fig. ^ we showed the percentage errors in obtained by two prescriptions for computing the derivatives. 
On the left we have used the implicit formula (|l^) , on the right we have simply computed a finite difference 
in O. We have also plotted the effects of adding more and more parameters into the problem. The lower 
curves invert only {Q, Q}. We then add more and more parameters, in the order in which they appear in s. 
This picture illustrates the uncertainties inherent in this formalism. It also shows how the errors may change 
dramatically by adding more and more parameters into the analysis. If one neglects massive neutrinos and 
uses a finite difference prescription to compute the curvature matrix then an error of 10% in the measurement 
of is obtained for w^ 1 — (15/i-ftT) 2 (deg) 2 and 9b — 0.5° (as in Q ||). By considering one massive neutrino 
this error goes up to 20% with the same prescription. If instead of using finite differences one computes 
dCi/dfl implicitly, the error goes up to 30%. 

5.2 The sample- variance/noise balance and undetermined parameters 

If one assumes the values of some of the parameters, so that fewer of them are allowed to vary, then sample 
variance naturally becomes less important. In this case, for example, it is no longer necessary to observe 
such a large number of fields to obtain a successful interferometric determination of SI. We may illustrate 
this point by considering the interferometer experiment discussed in the previous section, but this time 
fixing the values of all but four of the model parameters, and working only with the set {Q, SI, Sl^h 2 , h}. As 
shown in Figure ^, with less than 10 fields one could now obtain a better than 5% determination of SI for 
w^ 1 » (10/j/-r) 2 (rad) -6 . Moreover, it is clear that for relatively high noise levels, one could concentrate on 
a single field, to overcome the instrumental noise, and still obtain a useful result; this is because the sample 
variance in a single field in this context is only 5 — 10%. Thus, even for w^ 1 (50/iif ) 2 (rad)~ 6 , and for a 
single field, one could now measure SI with an error of about 10 — 15%. 
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Figure 9: Percentage standard errors in f2 for an interferometer with 9 W = 4°, considering only the 4 
parameters {Q, ft, ftbh 2 , h}. Isolines are plotted from 5% to 100% in increments of 5%. 



6 Interferometers and other cosmological parameters 

Although not central to our discussion, it is interesting to explore the marginal variances for cosmological 
parameters other than f2, in particular for interferometers. This is because interferometers may be expected 
to provide the least noisy estimates of the power spectrum at high I. Hence the detection of subtle but 
valuable high I signals, such as the acoustic oscillation signature, are obvious goals for such experiments. 



In Fig. 10 we show the percentage errors in the parameters f^/i 2 (top) and h (bottom), for an interfer- 
ometer with a primary beam width of 4° , and for which we have allowed all the CDM parameters to vary 
(but assuming no massive neutrinos). Estimates of the two parameters illustrated in the figure are derived 
mainly from the relative positions and height of the Dopplcr peaks in the CMB power spectrum, and it 
is clear from the figure that a reasonable determination of these parameters is possible. Assuming a noise 
level of w _1 w (lO/iAT^rad)" 6 , the error in both f^/i 2 and h is about 35%, even assuming such a large 
parameter set. As discussed above, however, these errors can be greatly reduced if the values of some of the 
other parameters, such as n s or A, can be fixed by other means. 



7 Ideal scales away from the vicinity of sCDM 

So far we have concentrated our attention on the errors associated in estimating f2 (and other parameters), 
assuming the sCDM values for the model parameters. In this section, we conclude by making a brief 
investigation of how the errors in the determination of Q may be much larger if the universe is not described 
by the sCDM model. 

In particular, we will consider two non-standard CDM models which illustrate two ways in which the 
problem of measuring fl may be more difficult. The first is an open universe with f2 = 0.3, which has all 
other model parameters set to their sCDM values; we call this model oCDM. If the signal is of this type then 
the ideal scales for weighing the Universe will be at much larger t. Consequently more stringent bounds will 
be put on the resolution of single-dish experiments. Interferometers, however, will be little affected. 

Our second model has Q — 1, but instead we assume a scalar tilt of ns = 0.8. The remaining parameters 
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Figure 10: Percentage standard errors in fli,h 2 (top) and h (bottom) for an interferometer with 6 W =4°, 
allowing all parameters to vary (but assuming no massive neutrinos). Isolines are plotted from 20% to 100% 
in increments of 5%. 
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Figure 11: Percentage standard errors in f2 for a satellite experiment with w^ 1 = (30fiK) 2 (deg) 2 in a universe 
with all parameters equal to those of sCDM, but with f2 = 0.3 (top) and ns = 0.8 (bottom). Isolines are 
plotted from 5% to 100% (top) and from 10% to 100 % (bottom), in increments of 5%. 

all have their sCDM values, and we call this model tCDM. For this type of model the intensity of the signal 
at the Doppler peaks is much smaller, requiring larger sensitivities from both inteferometers and single-dish 
experiments. This obstacle for an accurate measurement of 51 occurs also for signals coming from theories 
with low Q, non-zero r (a gravitational wave component), or non-zero r (late reionization) . 

7.1 Satellite experiments 

In Fig. |lj we plot the percentage errors in f2 for oCDM (top) and tCDM (bottom), for a satellite experiment 
with w^ 1 = (30fiK) 2 (deg) 2 . It is clear from the figure that for oCDM we require a resolution of better than 
0.4° in order to determine Q to any reasonable accuracy. Overall in comparison with sCDM the error plot 
gets squashed towards lower 6b. This is what one would expect, since the Doppler peaks in such a model 
are shifted to significantly higher multipoles, and so a finer angular resolution is required to determine their 
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position. For sufficiently fine resolution and with this relatively low level of noise, we also see that all-sky 
coverage is now not a great advantage. This is because the sample variance on the CMB power spectrum at 
the position of the Doppler peaks is reduced by their shift to higher multipoles. In any case, however, for a 
resolution of 0.1°, one could hope to achieve an accuracy of better than 5%, with little dependence on the 
area of sky covered. 

For tCDM, by comparing the bottom panels of Figs || and [ll], we see that again the errors are much 
larger. Overall the effect of reducing the power in the Doppler peaks by considering, say, a theory with 
Tig < 1, is similar to increasing the noise parameter For this reason all-sky coverage is now far from 

ideal. A 10% determination of under these conditions would now require a resolution of about 9t ~ 0.2° 
with an ideal scanning area corresponding to L as 100°. 

7.2 Interferometer experiments 

We have also investigated the performance of interferometer experiments in these two non-standard cosmolo- 



gies, and the results are given in Fig. 12 



From the figure we see that interferometers suffer far less in oCDM than satellite experiments. In fact, 
there is very little change in the accuracy to which O can be measured. If anything the situation improves 
wherever interferometers are limited by their large sample variance. This is not surprising, since the uv- 
coverage is assumed constant up to £ 1000. This still easily includes the relevant Doppler peaks in this 
model. Therefore the noise factor in the sensitive sections of the spectrum does not increase in the same 
way as it does for a single dish experiment. As an example, for w^ 1 w (50/iX) 2 (rad) _6 , we can expect an 
accuracy of about 15-20% by observing 10-20 fields. 

In the tCDM model, interferometers are affected in the same way as single-dish experiments. Again one 
suffers from a weaker Doppler peak signal, effectivelly similar to a larger noise parameter. Nevertheless, for 
i=s (10/iif ) 2 (rad)~ 6 with 20 observed fields, the error in is still around 20%. 



8 Summary and concluding remarks 

In this paper we reused the linear approximation tools developed in |B| in a broader context. We considered 
a larger class of experiments, more evasive signals, and formalism implementations producing less spectacular 
results. We then set off looking for the ideal experiment for weighing the Universe. The exercise is clearly 
somewhat naive. Political forces are more likely to determine experimental design than scientific reasons. 
Nevertheless it is interesting to note from our plots that accurate measurements of 51 may be expected even 
from poorly designed experiments. 

In a rough summary we found three corners in experiment parameter space which, for different reasons, 
and using different strategies, provide estimates for Q with useful errorbars. Satellite experiments make use 
of their low noise and intermediate resolution to map the first peak (and maybe part of the second) very 
accurately. They rely on beating sample variance by means of a large sky coverage. Ground-based and 
balloon borne experiments may make up for their larger noises with improved resolutions. This results in a 
better mapping of the second and maybe third Doppler peak. Small-sky coverage is used as a strategy to 
beat noise. The resulting increase in the sample variance is expected to be less of a problem given the larger 
£ at which the relevant measured features occur. Both these experimental strategies suffer dramatically if 
the Doppler peaks happen to be at larger £ than in the standard model, such as is the case with open models. 
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Figure 12: Percentage standard errors in Vt for the interferometer experiment with 6 W =4°, for a model with 
all parameters equal to those of sCDM, but with Vl — 0.3 (top) and ns — 0.8 (bottom). Isolines are plotted 
from 10% to 100% (top) and 10% to 100% (bottom) in increments of 5%. 
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Interferometers constitute a third viable avenue. They may be expected to map the first peak very 
poorly, but they will do outstandingly well in mapping the secondary peaks, even when compared with 
high-resolution single-dish experiments. Again a small coverage area technique is employed so as to reduce 
the effects of noise. The sample variance is always large, but again one may hope that its practical effects 
are small due to the high I features which should be targeted. However, this is not always the case, and 
in fact interferometers are often limited by sample variance. For this reason experimental interferometric 
parameters exist, for which performance is improved if the relevant features are at higher as with open 
models. All this seems to indicate that one should explore more than one experimental avenue in order to 
avoid unpleasant surprises. 

We conclude with a couple of remarks on how to improve the linear analysis performed here. Clearly 
one should compute the likelihood directly from the Ce(si) function. This is feasible but computationally 
demanding. Typically the non-linear analysis tends to predict smaller errors. Also one should go beyond 
the uniform prior assumption used here. A variety of information on all of the Sj parameters exists, coming 
from all sorts of other sources. This should be incorporated in the priors, naturally reducing the errors of the 
final estimates. For instance we have found errors in f2„ of order 0.9. These clearly are not compatible with 
what we already know about the mass of the neutrinos, should there be a heavy neutrino. The same can be 
said for some of the errors found for fib, h, ns- The final analysis waiting to be done should therefore not 
only go beyond the linear approximation, but also proceed to set up priors incorporating the vast amount 
of information coupled to the parameter determination problem posed by the CMB. 

In closing we should mention that marginal variances may be perhaps somewhat misleading. We have 
used them simply because they allow for a nicer display of results. However errorbars are in fact ellipsoids in 
the parameter space {s^}. In the non-linear approximation they become rather irregular bubbles of maximum 
likelihood. Such bubbles cannot be visualized but they encode all the information one may extract from a 
given experiment. Their computation and storage may turn out to be computationally very demanding, but 
such is the task one should perform when data finally starts pouring in. 
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